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Abstract: We study the effects of spacetime noncommutativity on the nonequihbrium 
dynamics of particles in a thermal bath. We show that the noncommutative thermal bath 
does not suffer from any further IR/UV mixing problem in the sense that all the finite- 
temperature non-planar quantities are free from infrared singularities. We also point 
out that the combined effect of finite temperature and noncommutative geometry has 
a distinct effect on the nonequilibrium dynamics of particles propagating in a thermal 
bath: depending on the momentum of the mode of concern, noncommutative geometry 
may switch on or switch off their decay and thermalization. This could have signifi- 
cant impacts on the nonequilibrium phenomena in the early universe at which spacetime 
noncommutativity may be present. Our results suggest a re-examination of some of the 
important processes in the early universe such as reheating after inflation, baryogenesis 
and the freeze-out of superheavy dark matter candidates. 
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1. Introduction 



On the noncommutative (NC) spacetime, the spacetime coordinates do not commute with 
each other anymore but obey the commutation relation: 

[x^'^ x"] = iO'"' (1.1) 

where 6'^" is a constant antisymmetric matrix |jl]]. Noncommutative quantum field theory 
(NCQFT) can then be derived from its commutative counterpart with the usual product 
of fields replaced by the Moyal star product: 

= et^'*^^M^^^0(x)x(l/)|.=.. (1.2) 
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As NCQFT arises naturally from string theory p|, it has received lots of atten- 

tions and has been an active research topic during the last decade. NCQFT has many 
unique properties such as Lorentz violation nonlocality and modified causality P, 0. 
Another intriguing phenomenon associated with NCQFT is the existence of the infrared 
(IR) /ultraviolet (UV) singularities p, ^ [1^. This is a phenomenon which gives rise to 
various pathologies in the field theory. Despite the loss of Lorentz invariance and local- 
ity, it has been argued that both CPT and spin-statistics theorems still hold |Tl|, |T2|, |13[. 
However, it has been pointed out that the space-time noncommutative theory (i.e. 7^ 
for i = 1,2,3) may violate unitarity if the theory also suffers from IR/UV mixing 



rR M. Therefore, to avoid getting into trouble with unitarity, we will confine ourselves 



to the case with 6'°* = in our study. 

On the other hand, nonequilibrium phenomena play a crucial role in many impor- 
tant processes in the early universe. These include reheating after inflation, baryogenesis, 
freeze-out of dark matter candidates, electroweak and QCD phase transitions fl^, 0, |19[. 
A common treatment of nonequilibrium evolution is to implement the closed-time-path 
(CTP) formalism [20, 21, 22, 2^ which is a path-integral approach to a time evolved 



density matrix. The thermal bath degrees of freedom are integrated out to obtain the 
nonequilibrium effective action which forms the generating functional for all the correla- 
tion functions. This approach also leads to quantum Boltzmann equations which can be 
solved to give the time evolution of the distribution functions. The accomplishment of 
thermal equilibrium is determined by the asymptotic time behaviour of the equal-time 
two-point correlation function and the distribution function. 

A natural question to be asked would be: what if we consider finite temperature 
and spacetime noncommutativity at the same time? Since noncommutative geometry 
naturally introduces a new energy scale i?NC ~ 9^^^'^ in addition to the temperature 
scale T, noncommutativity could have an interesting impact on the time evolution of a 
nonequilibrium system. In particular, if spacetime noncommutativity is really present in 
the early universe, it would be important to understand how does it affect the relevant 
physics. Some behaviours of NCQFT at finite temperature have been investigated in 



[p4| , |25| , |26| , P7| , P8| , |30| , pl[| . However, none of these works have considered the decay 
and thermalization of particles propagating in a noncommutative thermal bath. In order 
to fill the gap, in this paper we raise and investigate the following question: how does 
spacetime noncommutativity affect the nonequilibrium dynamics? 

In order to address this question, we consider a simple model with scalar $ particles 
propagating in a thermal bath constituted by two other different scalars Xi and X2- The 
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question of how do the $ particles come from at the first place is irrelevant to our discus- 
sion. Immediately after the $ particles are created, they are not in equilibrium with the 
thermal bath. They can either decay into Xi + X2 or thermalize with them, or both, de- 
pending on their kinematical properties. A complete understanding of the nonequilibrium 
dynamics of $ would be to study the time evolution of their correlation and distribution 
functions. However, as a first step, we will only confine ourselves to their decay and ther- 
malization processes in this study. The decay rate and thermalization rate Ft of the 
$ particles are characterized by the imaginary part of their self-energy as well as their 
in-medium dispersion relation |^2[. Notice that in the expanding early universe, even if 
Ft 7^ and thereby a thermalization process is kinematically favored, it still does not 
guarantee an actual thermalization. To maintain thermal equilibrium with the thermal 
bath in the early universe, we require Ti > H where H is the Hubble expansion rate. If 
Ft = or Ft < H, then the $ particles will continue to be nonequilibrium; and whether 
they can have out-of-equilibrium decays depends on Fd 7^ or Fd = 0. 

This article is organized as follows. In Section 2, we will outline our model and sketch 
how to compute the different contributions to the imaginary part of the self-energy, which 
will be analyzed to reveal its properties. We will also compute the real part of the 
self-energy and hence obtain the noncommutative in-medium dispersion relation. The 
IR/UV mixing issue in our model will be analyzed as well. In Section 3, we will study 
the impacts of spacetime noncommutativity on the decay and thermalization processes of 
the $ particles propagating in the thermal bath. We will find that as a combined result 
of finite temperature and noncommutative geometry, the stability and thermalizability 
properties of the $ particles propagating in the thermal bath are altered in a momentum 
dependent manner. Finally, in Section 4, we will give some preliminary discussions on the 
possible applications of our results to the early universe. 



2. The Model 



We consider a theory of three interacting real scalar fields in noncommutative spacetime 
with the following Lagrangian density 



1 1 ^ 



i=l 



|xi*$*X2- |x2*$*Xi, (2.1: 



where Mb is the bare mass of $. We will assume that the mutual interaction between the 
fields Xi ) X2 ensures them to be in thermal equilibrium at a temperature T = 1/(3. The 
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commutative counterpart of this model has been previously studied in ^ with an 



analysis of the different processes in the thermal medium. Here we will follow the similar 



treatment and conventions as [0 



The relevant quantity is the self-energy of the field $ which we will obtain to one loop 
order 0{g^) in the Matsubara representation. The one-loop self-energy is given by 

UJjn 

where Um = 27rm//5, Un = 2TTn/f3 are the bosonic Matsubara frequencies and p x k = 
Pi9^^ kj. The factors | and |e*''^'^ represent the planar and non-planar contributions 
respectively. Obviously the noncommutative phase factor is nontrivial only if fc* is nonva- 
nishing in the direction where 6'^^ is non-zero . The Matsubara propagators G^^^ and G^^ 
are written in the following dispersive form 

G^,{uj„,,p} = [ dpo ^'^^':^ , (2.3) 

J po- lUJm 

n I , - , 7"^^ [a P2(go,P + ^) 

G^^[uJm + Vn,P + k) = / dqo : : — , (2.4) 

where the spectral densities for xi ^'^'^ X2 are 

Pi(Po,p) = 77— ['^(Po-^i) -<5(po + ^i)], uji = J]P + Mf, (2.5) 

p2{qo,P + k) = -^[6{qo-uj2) -S{qo + uj2)], UJ2 = J {p + kf + M|. (2.6) 

Z L02 

This representation allows us to carry out the sum over Matsubara frequencies ujm in a 
rather straightforward manner [Rl, 53. The resulting self-energy can be further written 
in the dispersive form 

EK,fc) = -l fdcJ^^M^^i), ,2.T) 



where ImSR(co', k) is the imaginary part of the retarded self-energy which is defined by 
the analytic continuation 

T.Yi{h,k) = T.{i'n = -iko- t,k). (2.8) 
2.1 Imaginary part of the self-energy 

The retarded self-energy Sr(co') has cuts along the real axis. The discontinuity across 
these cuts (defined by Sr(ci; + ie) — ^^(a; — ie)) gives the imaginary part of Sr(u;): 

DiscSR(cu) = -2iImSR(cj). (2.9) 
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Xi,P 




X2,P + k 



X2,P + k 



Fig 1: One-loop processes contributing to aQ'^^{uj,k),a^'^^{u,k,T) and a^'^^ {u, k,T). 
The inverse processes are not shown. 



It is then easy to obtain 



lm^n{u,k)=g^f 



where 



(27r)3 2uJi 2uj2 



■ (1 + ^1 + ^2) 

+ b{u) + cui - uJ2){nx - n2) 
+ 6{uj + 002 - uji){n2 - Hi) 



Hi 



n\uji 



n[uo] 



(2.10) 



Here the factor | gives the planar contribution ImSj^, while the factor ^e^^^ gives the 
nonplanar contribution ImE^^. We will use the superscripts "P" and "NP" denote the 
planar and non-planar contributions respectively. 

It is convenient to write ImSj^''^^ as a sum of several contributions of different physical 
origin, namely 



lmT.^{uj,k,T) = aQ{uj,k) + al{uj,k,T) + al{uj,k,T) , 
ImE^^ico, k, T) = a^^'ico, k) + k, T) + K T) . 



(2.11) 
(2.12) 



Here the quantities Gq'^^ {uj, k) are the zero-temperature contributions, while crf'^^(co', k, T) 



a^'^^ {uj, k,T) are the finite-temperature contributions. At one-loop, the processes that 



contribute to aQ'^^{u}, k), a^'^^{uj, k,T) are $ ^ xi + X2, while the processes that con- 



tribute to (T^'^^{uj, k,T) are Xi,2 ^ $ + X2,i- See Figure 1. 

In the commutative version of our current model has been studied and the 
corresponding contributions to the imaginary part of the self-energy are computed as 
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(To, (Ta and ab respectively. In fact, ImS^ is precisely | times the corresponding results in 
32| , which implies that (Jq{uj, k) = ^ ctq, o"^(cc;, fc, T) = | (Jq, crf{u, k,T) = ^ ab- We have: 



ao{uj,k) 
al{uj,k,T) 



9' 



327rfc 

,2 



sign(cu) e[ - (Ml + Ma)' ]{B - A) , 



9 



327rPk 



32TrPk 



signHe[Q2_(Mi + M2) 



signHe[(Mi-M2)2-g2 



In 
In 



,-f3B 



1 - e-f^^ 



'HA 



+ (Mi 



(Ml 



M, 



(2.13) 
, (2.14) 

, (2.15) 



where A 



B 



(2.16) 



and 



\u;\[Q^ + {Ml - M|) ] ± fcy/jg^ + (Mf - M|) - Mf 

2g2 



(2.17) 



In 0"^ and , l3k <^ 1 corresponds to the high temperature limit while f3k ^ 1 corresponds 
to the low temperature limit. Both and increases with decreasing jSk. For instance, 
cr^^ when pk oo. On the other hand, both and af approach to a finite value 
when jSk —>■ 0. 

As for the non-planar parts, we have to include the factor e*^^'^ in the integral. To 
evaluate the integral, we note that since the vector 6^^kj is perpendicular to fc*, it is 
convenient to adopt a spherical coordinate system with the polar axis pointing in the fc* 
direction and the azimuthal angle measured from an axis defined by 6'^^kj. Denoting 
the polar angle by we obtain 



p X k = p9 k sm'd cos ( 



(2.18) 



where 6 is the magnitude of the vector O^^kj and A;* is the unit vector of fc*. For example, 
if /c* = (0,0, fc), then 6 = Ay^(^^^0^~+~(^^^- The calculations then proceed in a similar 



fashion as has been performed in the appendix of After lengthy but straightforward 
calculations, we obtain the non-planar contributions to the imaginary part of the self- 
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energy 

'Jo'' = ^ sign(a;) Q[Q^-{M^ + M2f]^ sin {il kQ 6) , (2.19) 

' = iC; ^is^(^) e[ - (M, + M,)^ ] '"t^''"' sin (fi . /Y^y^) 

^ V V V / / 

+ (Mi^M2), (2.20) 

©(-<5^) . sinh {n Pk n) 



NP y • t \ 



n=l 



+ e(g^)e[(M,-M,)^-Q^] E '/T, sin ^ 

+ (Mi^M2), (2.21) 

where Q = (S — A)/(2A;) in the above expressions^. We note that \o^^\ < \af\ for 
r — 0,a, b. 

Prom the above expressions, it is obvious that a^^ — > (Jq when kQ9 <IC 1. In fact, to 
the leading order, Q ~ Mq where Mq is the renormahzed mass of $ (to be determined in 
the next subsection). Thus, kQ9 <^ 1 corresponds to the case where 9 <^ (I/Mq) (l/k). 
Notice that V9 = Anc represents the characteristic length scale smaller than which the 
effect of noncommutative geometry becomes significant. Moreover, the Compton and de- 
Broglie wavelengths associated with $ can be identified as Ac ~ I/Mq and Ads ^ l/k 
respectively. As a result, we conclude that Cq when Anc ^ V^c AdB- This is 

always true when Anc is much smaller than the smaller of Ac and AdB- In the relativistic 
case, AdB ^ Ac; while in the non- relativistic case, Ac <^ AdB- In any case, if Anc is really 
much smaller than the smaller of Ac and AdB, then the resolution due to either Ac or Ada 
is not high enough to see the effect of noncommutativity. The system behaves as if it 
were commutative, and hence cr^^ — > a^. On the other hand, a^^ when kQ 9 » 1, 
which corresponds to the case Anc ^ "s/AcAdB- This is always true when Anc is larger 
than the larger of Ac and Ada- In this case, the system is completely noncommutative. 

For a^^ and af^, the significance of noncommutativite effect depends on the ratio 
Q9/P ^ MqT9. When Q9/I3 <^ 1, noncommutativity is neghgible. This corresponds 
to the case when Anc <S At where At ~ 1/ \/MqT can be identified as the thermal de- 



^Notice that if {Q9/I3f - < 0, we replace sin(...) by sinh(...) and ^J{Q 9/P)'^-n^ by 
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Broglie wavelength associated with $ when it is propagating in the thermal bath with 
temperature T. In the thermal medium, it is the thermal de-Broglie wavelength that 
plays the role of the characteristic resolution acquired by When Anc ^ '^t, the 
system behaves as if it were commutative, and hence cr^^ — > cr^^. On the other hand, the 
effect of noncommutative geometry is significant when QQj^i^ 1, which corresponds to 
Anc ^ At- In this case, both a^^ and a^^ are suppressed. 

In general, as Q increases, both and decrease. This can be attributed to 
the reduction of degrees of freedom due to spacetime noncommutativity. In particular, 
as ^ — > oo, cr^^ ^ 0. In this case, the reduction of degrees of freedom due to spacetime 
noncommutativity is maximal. This can also be understood from the extremely rapid 
oscillations of the phase factor in the integrands of and a^^ . 

Furthermore, similar to the commutative case, i5k <^ 1 corresponds to the high tem- 
perature limit while ^ 1 corresponds to the low temperature limit. This is true 
regardless of the relative significance of the noncommutative geometry. Both a^^ and 
(7^^ increase with decreasing fik. For instance, a^J — > when /3A; — > oo. On the other 
hand, both a^^ and a^^ approach to a finite value when pk — > 0. 



2.2 Dispersion relation 



The real part of the self-energy is given by 

ReEK(.„, k) = p. + M? + Uf + mU + *^ ^ ^ ■ 

Again, the factors | and \ e^^^^ represent the planar and non-planar contributions respec- 
tively. To facilitate the calculation, we introduce the Schwinger parameters 

1 



1 



POO 

-- / (iaie-"i(^"'+^?+'^-), (2.22) 
Jo 

poo _ 
{p + k)^ + Mi + {uJm + Un)^ Jo 

By completing squares, the p integrals now becomes Gaussian and can be readily 
evaluated to give 

ReE^'^Pfi/ k)^-^ r ^ e-4i4+i(^i+^l)]-^ 

• r (ia:e-3"'^-^'-5"(^?-^2)-^?(' !!(l-a;), ^V(2.24) 
J-i V 2 47rQ; / 
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where a = ai + 0^2, x = (ai — 02) /a, k'^ = k'^ + z/^, and 



k"^ 

c 

E 



2iTimz+iiTm'^ 1 



(2.25) 



is the Jacobi theta function |3^. Similar to §], we have multiphed the above integrands 
by exp (— 1/a A^) in order to regulate the small a divergence such that 



1 



for Re 



-1 + ^, for ReSNP ' 



(2.26) 



where k = kO. 



The leading contribution of the integral ( p.24|) comes from the region a ~ 0. After 
performing the above integrations and upon the analytic continuation i/„ —iuj — e, we 
obtain the retarded self-energy Re ^^(ct;, k) = Iq + Iq^ + 1^ + 1^^ with 



tP, np 
'0 



rP, NP 



9' 



327r2 
2 



9 



InfcL 



9'Q' 



00 _ 

E A 

Z — / 

m=l 



for cL < 1, 

- M| 



(2.27) 
(2.28) 



1 + 



+ 1 



Ml 



Ki{2MiL') 



Ml - M| \ L' 



i^i ( 2 M2 L' ; 



(2.29) 



where Q'^ = uj'^ — k"^ = —k\. Kq and Ki are the modified Bessel functions of the second 
kind. Here 

, (g2 _ - M^f) ((Ml + M,f - Q^) 



and 



'12 



(2.30) 



(2.31) 



In fact, Iq'^^ and J^'^p arise from the m = and m 7^ terms in the Jacobi theta 
function respectively. In the zero temperature limit T — > 0, the only non-vanishing term 
in the Jacobi theta function comes from m = in which case i){z, r) ^ 1. Notice that c 
is always positive-define if (Mi - Ma)^ < < (Mi + Ma)^ For > (Mi + M2Y, the 
$ particle can decay into xi and X2- For Q'^ < (Mi — M2)^, Xi (or X2) can decay into $ 



- 9 - 



and X2 (or xi)- In both cases, c becomes imaginary and a non-zero imaginary part will 

7-P, NP 

appear m Jq 

The bare mass Mb of $ receives renormalization from both of Iq and Iq^, and so 
the (zero-temperature) renormalized mass for $ is defined as 

Mq = MI + Iq\q2^m§ + Io^\Q^=Mg- (2.32) 
To the order 0{g^), the dispersion relation is then given by 



2- /t^ + M^+jP+jNP, (2.33) 



where 

T'' ^ IX.^M^^ and = 7^nQ^=M„^ (2-34) 

represent the finite-temperature corrections to the dispersion relation. 

It is instructive to examine the behaviour of the finite temperature quantities 
in various limits. In the low temperature limit with T -C Mi, M2, we have M12L' ^ 1. 
Using Ki{x) \f^^'^ for X 1, it is obvious that I^^^^ are exponentially suppressed 
as long as T <^ Mi, M2, regardless of the significance of the noncommutativity. This is 
consistent with the fact that 7^' arise from the m 7^ terms in the Jacobi theta function 
and can survive only at finite temperature. Similarly, T^^ is exponentially suppressed in 
the large Q limit with 1/kd <^ Mi, M2, regardless of the magnitude of the temperature 
T. On the other hand, in the high temperature limit with Mi, M2 ^ T, l/kO, we have 
Mi^2L' -C 1. Using Ki{x) ^ - for x <C 1, we obtain 



rP,NP _ 9^Q^ 



96 {(3u;y 



^ ^ Mf-M| \ 1 ^ / M!-MI\ 1 



(2.35) 



It is remarkable to notice that the above expression for /P-'^^'p are independent of m and 
the ultraviolet cut-off A. Also, it is vahd irrespective to the relative magnitude between T 
and l/kO. Most interestingly, the non- planar contribution I^^ is completely independent 
of 9. This means that in the hmit Mi, M2 ^ T, l/kO, the dispersion relation receives 
vanishing finite-temperature corrections from the spacetime noncommutativity. 

Finally, let us comment on the IR/UV mixing cfi'ect of the noncommutative field 
theory at finite temperature. Needless to say, since Iq^ is a zero-temperature non-planar 
quantity, it does suffer from the usual IR/UV mixing. Since the IR singularities are a 
refiection of the fact that the field theory is UV divergent, the key to resolve the IR 
singularities lies at a proper UV finite completion of the noncommutative field theory 
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p7| . With a choice of UV completion, the IR singularities will get smoothen out. For 



example, a natural choice is to embed the noncommutative field theory as a low energy 



field theory of open string theory in background i?-fields. It was shown explicitly in ^7 
that in doing so, the IR pathologies of noncommutative field theory are resolved. In 
particular, in the deep IR, the theory fiows continuously to the commutative field theory 
and the normal Wilsonian behaviour is restored. Therefore, with this understanding, the 
zero-temperature IR/UV singularities are harmless. 

It is clear that our noncommutative thermal bath does not suffer from any further 
IR/UV mixing problem, in the sense that all the finite-temperature non-planar quantities 
are healthy and are absent from any infrared singularities if we take the limit k ^ after 
taking the limit A — > oo. For instance, let us look at the quantity L' in J^^ (which is of 
purely finite-temperature nature) as displayed in ( p.29|) . It is obvious that L'^ 
when A — > CX3 and — > 0. Unless T oo, I^^ is manifestly finite. Physically, it is 
the finite temperature T acquired by the thermal bath that rescues the system from any 
further IR/UV mixing problem: the finite temperature T acts as an effective ultraviolet 
cut-off for the system once we have taken the limit A — > cxd. 



3. Effects of Noncommutativity on Nonequilibrium Dynamics 

In previous section, we have computed the imaginary part of the self-energy as well as 
the dispersion relation associated with the $ particles. Now, we can proceed to study 
the nonequilibrium dynamics of $ when it propagates in the noncommutative spacetime. 
In particular, we will focus on how does spacetime noncommutativity affect the possible 
decay and thermalization processes of $ in the thermal bath. Without loss of generality, 
we will assume that Mi > M2. 



3.1 Decay and thermalization rates 



Let us write the imaginary part of the self-energy as ImE^ = ctd + ctt where ud = 
ctq + 0"^ + 0"^^ + cr^^ and o"t = o-^ + a^^. These quantities are regulated by kinematical 
constraints and take the forms: 

= e[Q2_(Mi + M2)'] AD(Q',fc2), (3.1) 
ax = e[ (Ml - M^f - Q' ] AT(Q^ + e( -Q^ ) A'^{Q^ k'), (3.2) 

where Ad{Q\ k'^) can be read off from ( p7[3|) , (g3§, {^J§, (g), while At{Q\ and 
Aij.{Q'^ , k"^) can be read off from ( |2.15| ), ( |2.21| ) respectively. To compute the decay and 
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thermalization rates |3^, we need to put the $ field on-shell, i.e. setting 

Q2 = m2+xP+JNP. (3.3) 

Then, to the order 0{g'^), we obtain 

Ad(M2,A;2^ 



D 



and 



At(M2,P) 



Q[Q' -{M1 + M2Y] (3.4) 



Q[{M,-M2Y-Q'l (3.5) 



where uq = ^Jk"^ + Mq . Note that up to the order 0{g'^), it is sufficient to set = Mq 
in Ad and At; however, one has to use the full expression (^.3]) inside the B functions. 



Physically, ctq + o"^^ represents the zero-temperature planar + non-planar contri- 
butions to the decay rate, while + a^^ represents the finite-temperature planar + 
non-planar contributions to the decay rate. Similarly, cr^ + a^^ represents the finite- 
temperature planar + non-planar contributions to the thermalization rate. Notice that 
af + a^^ is a purely finite-temperature effect and there is no zero-temperature counterpart 
of it, namely + af^ — > as T ^ 0. 

The rates T^lk) and TT{k) depend on k, the magnitude of the spatial momentum. 
If TY){k) = for a certain k, then this mode of the $ particle is stable in the thermal 
medium. Conversely, if T]:,{k) 7^ for a certain k, then this mode would be able to decay 
into Xi + X2- Kinematically, the feasibility of decay is regulated by the Heaviside function 

e[Q^ -{Mi + M2f]. (3.6) 

On the other hand. Ft (A;) 7^ for a certain k implies that this mode of $ acquires a 
relaxation or thermalization time scale beyond which it approaches thermal equilibrium 
with the bath constituted by xi and X2- This happens through the decay of xi,2 into 
X2,i and $, and their recombination, namely xi,2 ^ X2,i + ^ [^H, 0- As the modes 
carrying momentum number k propagate through the thermal bath, they will be screened 
or dressed by the excitations in the medium and will propagate as quasi-particles. In 
fact, TT^{k) characterizes the "decay rate" of the quasi-particles associated with $ in the 



medium, and this is precisely the relaxation or thermalization rate ||3^, q^. Of course, 
if rT(A;) = 0, then it takes an inffnitely long time for these modes of $ to approach 
thermal equilibrium with the bath, which simply means that these modes can never 
thermalize with the bath. Kinematically, the feasibility of thermalization is regulated by 
the Heaviside function 

e[{M,- M2Y -Q^]. (3.7) 
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3.2 Momentum dependent alternation of stability and thermalizability 



What we would like to study is the following. How does the conventional understanding 
of the dynamics of $ in the thermal bath change in the presence of noncommutative 
geometry? 

Since the decay and thermalization processes depend crucially on the relative size of 
compared to (Mi — M2Y and (Mi + M2)^, it is useful to consider the following three 
regions of Q^: 



7^l 
7^2 
7^3 



Q'^ < (Ml - Ma)^ where we have Ft 7^ and To = 0, 

(Ml - Ma)^ < < (Ml + M2)^ where we have Ft = Td = 0, 

(Ml + M2f < Q^ where we have Ft = and To 7^ 0. (3.8) 



If we define Qq as the commutative counterpart of Q^, then it would be given by Qq = 
Mq + 2X^. When one turns on the noncommutativity parameter, Qq changes to and 
could cross over from one region to another, corresponding to turning on or turning off 
the decay and/or the thermalization processes as one takes into account of the effects 
of noncommutativity. What can actually happen depends on the "jump" — = 
jNP _ jp gi^pg |jNP| < |jP|^ always have _ g2 ^ q if > and 

— Qq > if < 0, regardless of the signs of X^^ . Therefore the direction of the 
"jump" depends solely on the sign of X^. By setting A — > 00 and Q^ = Mq in ( f^.29| ), it 
is easy to obtain that 

X^>0 when MS<^-^§0f, (3.9) 

X^ < when M^ - M^ < M^ (3.10) 

for any temperature T and momentum k. And for Mq in the intermediate region, one 
obtains 

XP>0 forT<To,l , (Mf-MlY (3.11) 

} when ^-j. < Mq2 < Mf - M|. ^ ^ 

XP<0 forr>To,J M2 + M| - °- 1 ^ 

Here Tq is a temperature that is determined by solving X^(To) = 0. Note that Tq depends 
not just on the mass Mq of the propagating particles $ and the properties of the thermal 
bath (i.e. Mi,M2), but also on the momentum k of the mode of $ in consideration. 
Therefore, as we remarked above, the impact of spacetime noncommutativity on the 
decay and thermalization processes is going to be momentum dependent. We remark 
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that due to the rather comphcatcd form for X^, an analytic expression for Tq is not 
available. Fortunately, this is actually not needed in our analysis below which is aimed 
at explaining the general physical features. In a concrete phenomenological study, a more 
detailed knowledge of Tq maybe needed and this can always be obtained numerically. 

As we mentioned above, the decay and thcrmalization processes may be turned on 
or off when one takes into account of the effects of noncommutativity. This happens 
whenever crosses over different regions as noncommutativity is turned on. As a result, 
we obtain the following six interesting cases which may happen in principle. 

1. Qlin region TZn Ql < (Mi - Ma)^ 

In the commutative picture, we have Fd = and Ft 7^ which imply that $ is stable 
and cannot decay into Xi + X2, but it can thermalize with the bath. Once we switch to 
the noncommutative bath, we may have modes such that 

• Case A: (Mi - Ma)^ < < (Mi + Ma)^ =^ Fd = Ft = 0. This corresponds 
to the situation with < 0. These modes of $ are still stable against decay into 
Xi + X2; but they cannot thermalize with the noncommutative bath anymore. In 
this case , we see that noncommutative geometry suppresses thermalization for these 
modes. 

• Case B: (Mi + Ma)^ < Q'^ =^ Fd 7^ and Ft = 0. This also corresponds to 
the situation with X^' < 0. These modes of $ which are originally stable in the 
commutative bath, can now decay into xi + X2 in the noncommutative bath, but 
they cannot thermalize with it anymore. In this case, noncommutative geometry 
induces decay but suppresses thermalization for these modes. 

2. region 7^2: (Mi - M2)^ < Qc < (^1 + ^2)^ 

In the commutative picture, we have Fd = Ft = which implies that $ can neither 
decay into xi + X2 nor thermalize with the bath. Once we switch to the noncommutative 
bath, we may have modes such that 

• Case C: < < (Mi - Mif =^ Fd = and Ft ^ 0. This corresponds to the 
situation with X^ > 0. These modes of $ still cannot decay into Xi + X2, but they 
can now thermalize with the noncommutative bath. In this case, noncommutative 
geometry induces thermalization for these modes. 
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• Case D: (Mi + Ms)^ < =^ To ^ and Ft = 0. This corresponds to the 
situation with < 0. These modes of $ can now decay into Xi + X2 in the noncom- 
mutative bath, but still cannot thermalize with it. In this case, noncommutative 
geometry induces decay for these modes. 



3. Qlin region TZsi (Mi + Afa)^ < Qc 

In the commutative picture, we have Fd 7^ and Ft = which imply $ is unstable 
and can decay into Xi + X2, but it cannot thermalize with the bath. Once we switch to 
the noncommutative bath, we may have modes such that 



• Case E: < < (Mi - Ma)^ ^ Fd = and Ft ^ 0. This corresponds 
to the situation with > 0. These modes of $ cannot decay into xi + X2 in the 
noncommutative bath anymore, but they can now thermalize with it. In this case, 
noncommutative geometry suppresses decay but induces thermalization for these 
modes. 

• Case F: (Mi -Ms)^ < < [Mi + M2f =^ Fd = Ft = 0. This also corresponds 
to the situation with X^ > 0. These modes of $ cannot decay into Xi + X2 in the 
noncommutative bath anymore, and they still cannot thermalize with it. In this 
case, noncommutative geometry suppresses decay for these modes. 



For a given thermal bath and a propagating particle, the set of masses Mq, Mi, M2 
and the coupling g are fixed. To decide which of the above six cases can occur, one 
need to examine carefully how changes when noncommutativity is turned on. As an 
illustration, we present the detailed analysis of cases A and B in the appendix A. The 
result is that cases A and B can occur only if the masses of the thermal bath satisfy 



-^>4 + ^. (3.12) 



and if Mq falls in the range of 

(Mf - M|)^ 
Mf + M| 



< Mq^ < -(Ml -Ms)^ (3.13) 



When these conditions are satisfied, then depending on the momentum k of the mode of 
concern, either case A or case B or both cases can occur. See ([A.17|) -( |A.23D . 
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Similar analysis can be performed for the other four cases. The conclusion is similar. 
For cases E and F to occur, the necessary conditions are that the masses of thermal bath 
has to satisfy 

M1/M2 > 11 (3.14) 

and Mq has be in the range 

^(Mi + Ms)^ < < Ml - Ml. (3.15) 



For cases C and D, no condition is needed on M1/M2. But ( |3.15|) is needed for the case 



C and ( |A.6| ) is needed for the case D to occur. Similar conditions as ([A.17|) - ([A.23|) can 



be written down for these four cases. We will skip them here. 

All in all, we conclude that all of the six cases A-F listed above can occur in general. 
Depending on the masses, temperature of the thermal bath and the magnitude k of the 
spatial momentum of the mode of concern, the decay and thermalization rates of $ can be 
either induced or suppressed when one takes into account of noncommutativity properly. 
In general, there is a window of momentum modes whose stability and/or thermalizability 
properties are altered. The window depends on the temperature of the thermal bath. 
Typically, the affected modes are in the high (respectively low) k regime when T is low 
(respectively high). 



4. Conclusions 

In this article, we studied the nonequilibrium dynamics of a scalar $ propagating through 
a noncommutative thermal bath. We showed that noncommutative geometry has a dis- 
tinct impact on the nonequilibrium dynamics of particles propagating in a thermal bath 
by providing a momentum dependent enhancement or suppression of their decay or ther- 
malization processes. This is a combined effect of finite temperature and noncommutative 
geometry. Also, we pointed out that the finite temperature T of the thermal bath can 
play the role as an effective ultraviolet cut-off which rescues all the finite-temperature 
non-planar quantities from any further IR/UV mixing problem. 

Although our thermal bath is represented in a specific way by the two scalar particles 
Xi and X21 our analysis is general and our conclusions should apply generally for any 
bosons in contact with a noncommutative thermal bath. 

A particularly interesting arena of application of our results is the early universe at 
which spacetime noncommutativity may be present. Our results call for a re-examination 
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of some of the important processes in the early universe such as reheating after inflation, 
baryogenesis and the freeze-out of superheavy dark matter candidates, which are generally 
believed to occur at very high energy scales where spacetime noncommutativity could be 
significant enough. It will be very interesting to re-examine these processes in details and 
see how the presence of noncommutative geometry could affect the conventional picture. 
This could also provide another channel for probing the presence of noncommutative 
geometry in the early universe, in addition to the possible cosmological imprint left by 
inflation, see for example PSI, BUI EO]. 
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A. Analysis of cases A and B 

Modes for which the cases A and B can occur have to satisfy the following conditions: 



Ml + 2JP < (Ml - M2)^ (A.l) 

|JP| <0.1M2, (A.2) 

(Ml - M2f < M2 + JP + < (Ml + M2)^ (A.3) 

(Ml + M2f < M2 + JP + . (A.4) 

Here the first condition specifies that one is initially in the region TZi. The second condition 
( |A.2| ) is a condition on the size of the thermal correction to Mq. We have imposed a 
conservative 10% correction so as to guarantee that we are in the validity regime of the 
perturbation theory. In principle, one could have taken a different value for the RHS of 
(|A.2|) , e.g. 0.2Mq. This will not modify the analysis below and the physical effects we are 
going to point out are generic. The condition ( |A.3|) is for case A and the condition ( |A.4| ) 
is for case B. 

To analyze these conditions, we first note that cases A and B can occur only if < 0. 
This means that it is necessary to require Mg > (Mf — M|)^/(Mf + M|). On the other 
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hand, the conditions (|A.1D and ( |A.2| ) imply that 



-[Ml - (Ml - M^f] < -JP < O.IMqI 



(A.5) 



This is a non-empty condition for only if Mq < |(Mi — M2Y. This condition, together 
with the previous condition on Mq , gives 

{ Ml - Ml 
Ml + M| ~ ~ 4 



2\2 



< Mq < -(Ml - M2 



(A.6) 



This specifies a non-empty range for M^ only if (Mf - M'^f/{Ml + M|) < |(Mi - M2) 
i.e. 

/\/f, 

(A.7) 



Ml ^ 

M2 



The condition ( [A.7| ) on the ratio M1/M2 and the condition ( |A.6 ) on Mq are common to 
cases A and B and are the necessary conditions for the cases A and B to occur. Now 
assume (|A.7] ) and ( |A.6|) are satisfied, we are guaranteed that ( |A.5| ) specifies a non-empty 
range for — X^. The question is whether and for what configuration of modes will ( |A.5|) 
be satisfied. To address this, we need the form of as a function of k and T, which we 
collect here (also for X'^^): 



X^ 



T 



NP 



m=l 

^2 A/f 2 °° 

5' ^^^0 



2 /\^2\//\^2 



m=l 



m2Mi2 
_m2 Mf 



/;,(m Xi 



h{myi) + 



m2M| 
m?Ml 



/;,(m X2) 



h{my2) 



Here c± = 1 ± (Mf - MD/M^ 



Mi/T, Vi = yiTr^Tw M,, ^ = 1, 2 



(A.8) 
(A.9) 

(A.IO) 



and 



h{x) := xKi{x). (A.ll) 

To analyze which set of modes satisfy the condition ( |A.5| ), we note that — X^ = at T = 
and -XP ~ d X T2 for large T > Mi, M2, where d ^ g^M^ /{96lu^)[c+/MI + c^/M^ > 0. 
Therefore, for any given k, there is a temperature range 

T^{k) <T <T2{k) (A.12) 

such that ([A.5| ) is satisfied. Here the temperatures Ti 2(A;) are determined by solving 



XP(A;,Ti(A;)) = -[M2-(Mi-M2 



and XP(A;,T2(A;)) = O.IM^. 



(A.13) 
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-JP(A:), k fixed 



0.1M2 



|(M2 - (Ml - M2)2) 




Fig 2: Solving the condition ( [A .51) . 



So far tlie obtained conditions ( A.7 ), (|A.6|) and ( |A.12|) are common to botli case A 
and case B. Let us now proceed to analyze tlie conditions ( |A.3| ) and ( [A.4|) . Define the 
temperatures Ti 2 by 



-(jP+X^P)(A;,Ti(A;)) = Ml-iM^+M^f and -(XP+J^P)(fc, r2(A;)) = Mo'-(Ml-M2)^ 

(A. 14) 

then for a mode with momentum A; in a thermal bath with temperature T, 



case A occurs if Ti{k) < T < T2{k) , 
case B occurs if T < Ti{k). 



(A.15) 
(A.16) 



Now since -(JP + J^P) < -2TP, it follows that T2{k) > Ti{k) for all k. This is important 



since it implies that the condition (|A.12|) has always a non-trivial intersection with the 
conditions ([A.15| ) or ( [A.16| ). For each mode with a fixed momentum k, depending on the 
relative sizes of Ti{k),T2{k),Ti{k),T2{k), either case A or case B or both cases can occur. 
Specifically we obtain the following result: 



When f2{k) < T2{k): 



if Ti(fc) < Ti(fc), then 
if ri(A;) > Ti{k), then 



case A occurs for Ti{k) < T < T2{k); (A. 17) 

case A occurs for Ti{k) < T < T2{k); (A. 18) 

case B occurs for Ti{k) <T <fi{k). (A. 19) 



- 19 - 



When f2{k) > T2{k): 



if fi(fc) > T2{k), then case B occurs for Ti{k) < T < T2(A;); (A.20) 

if Ti{k) < fi{k) < T2ik), then case A occurs for fi{k) <T < T2{k); (A.21) 

case B occurs for Ti{k) < T < fi{k); (A.22) 

if fi(fc) < Ti{k), then case A occurs for Ti{k) <T <T2{k). (A.23) 

This means that given the temperature T of the thermal bath, modes with different 
(spatial) momentum k can behave quite differently in regard to the thermalization or sta- 
bility properties. In particular only a certain window of momentum modes as specified by 
the inequalities in ( |A.17| )- ( |A.23D above will be affected for a given T. Since the character- 
istic temperatures Ti^2{k),Ti^2{k) are monotonic decreasing function of k, it means for low 
temperature T, the affected window is in the higher k regime; while for high temperature 
T, the affected window is in the lower k regime. This momentum dependent alternation 
of the stability or thermalizability is novel and is a combined effect of finite temperature 
and noncommutative geometry. 
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